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• In this paper we introduce a diagnostic for measuring the quantum-classical difference for open 

' quantum systems, which is the normalized size of the quantum terms in the Master equation for 

, Wigner function evolution. For a driven DufHng oscillator, this measure shows remarkably precise 

scaling over long time-scales with the parameter <^o = /D. We also see that, independent of C,o 
the dynamics follows a similar pattern. For small all of our curves collapses to essentially a single 
curve when scaled by the maximum value of the quantum-classical difference. In both limits of large 
and small we see a saturation effect in the size of the quantum-classical difference; that is, the 
instantaneous difference between quantum and classical evolutions cannot be either too small or too 
CN '. large. 

PACS numbers: 05.45.Mt,03.65.Sq 
^ ■ 

The quantum-classical transition for open quantum systems is important to understand for fundamental and 
C ' practical reasons, including the design of quantum computers. For quantum systems with classical analogues permit 
chaos this transition is particularly interesting because these systems often display unusual quantum effects, such as 
rapid entanglement generation [l[ and hyper-sensitivity to perturbation Q , and it would be useful to understand when 
such effects appear. It is an inherently multi-parameter transition depending on the relative size of Ti compared to the 
characteristic action and on the strength of the interaction between the system and its environment as measured by 
some parameter D. The chaotic behavior of the classical limit is also crucial, as seen in their decoherence dynamics 
' argued to be related to the Lyapunov exponents of the system. An efficient approach that allows for broader conclusions 
lO ■ is to use Q scaling in composite parameters, which indicates @ that the quantum-classical difference as measured by 
' some quantity QCdi^T-, D, A) is a function QC'^{C) of a single composite parameter ( = h'^D^X' . 

, A good measure QCd is intuitive, global, and easy to compute. A KuUback-Liebler-like distance between classically 
' propagated and quantally propagated distributions picks out scaling properties cleanly Q when computed for a 2- 
, phase-space dimension map. Overlaps functions are however, difficult to compute in greater dimensions, or for flows. 
Studying the quantal dynamics entropy iS'2 = \n.{P), where P = Tr(p^) also yielded useful insights. However, the 
scaling did not last very long in that instance, and the direct relationship to quantum-classical difference is not clear; 
the use of the entropy has also been otherwise critiqued We report here on a measure that is (a) time-dependent, 
(b) explicitly measures the difference between quantum and classical evolution, and (c) is calculated from the evolution 
^ ' of a single distribution. This measure shows remarkable scaling, in parameter space and over long time-scales, and 
we consequently uncover interesting insights into unexpected behavior in the dynamics. Specifically, in the classical 
limit, the relative size of the quantum and classical terms saturates, so that quantum-classical differences continue to 
be propagated, rather than decreasing with time as might be naively expected. Conversely, at the near-quantal limit 
of small D and large Ti the quantum terms remain comparable to the classical instead of dominating. Taken together, 
this quantifies the smoothness of the quantum-classical transition for open systems compared to closed systems. 

We start by considering the Master equation for a Wigner function pw, or quantum quasi-probability, evolved 
under Hamiltonian flow with potential V{q) while coupled to an external environment Q: 

, , d^-+'V{q) d^-+'pw , 

= {H, pw} + 22"(2n-Hl)! Sg2„+1 dp2n+l + D^PPW (1) 
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The Poisson bracket Lc generates the classical evolution for pw, the quantal Ti terms are denoted by Lq and the 
environmental coupling T is modeled by a diffusive term with coefficient D. The computational results presented 
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below use coupling only to the momentum variables. For analytical simplicity, we assume coupling to all phase-space 
variables, justified since the dynamical chaos mixes the various phase-space directions. 
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FIG. 1: Quantum-classical difference as measured by {G{t)}, which is G{t) averaged over one driving period, for various initial 
parameters. The y-axis is logarithmic. Notice the remarkable scaling with the composite parameter ~ /D. For the initial 
Gaussian to be well-localized in the chaotic region, < 0.05. Hence for our calculations, we set o", = 0.05 and determine ffp 
from the constraint, aqap = h/2, imposed by minimum-uncertainty condition 

The measure we propose in this paper is a normalized average of the square of Lq (see Eq.( [1} for the precise 
definition). The square before the averaging is necessary to compensate for the arbitrary sign of Lq in different 
regions of phase-space. For example, in the driven Duffing problem that we study below and which is given by 
H = /2m — Bx^ -\- [C /2)x^ + Axcos{ijjt) the quantity Tr[Lq] is identically 0. The normalization makes the measure 
dimensionless, yielding 

G{t) = (2) 

Physically G is the relative size of the quantum part of the evolution for the Wigncr function, and it can exceed 
unity through partial cancelation between Lq and the Poisson bracket in the evolution equation for p^. Note that 
G does not measure the classicality of a state, but rather the classicality of its evolution and is analogous to the 
time-derivative of previous-used measures. 

In Fig.[TJ we plot G{t) for various combinations of h and D for a Duffing oscillator with the parameters (m = 1,5 = 
10, C = 1,^ = 1, and uj = 5.35). These results are calculated starting with a minimum uncertainty Gaussian Wiener 
function that is well localized in the chaotic region and propagating it using Eq. ([T]) as has been previously doneQ 0]- 
The figure shows remarkable scaling behavior in the measure G(t). That is, the dynamics of the quantum-classical 
difference are seen to depend only on the composite parameter Co = ^ 1^ over a factor of 800 in its value. Although 
such behavior has been seen earlier for short times with the entropy 0] , in this case the scaling dynamics lasts for 
the full time-scale monitored by us. 

The absolute size of G grows with a-nd is much greater for Co — 40 compared to Co — 0.05, which is reassuringly 
physically intuitive. Specifically, the classical limit is the regime when Co ^ 0.2, where G is always so small that the 
differences between the evolution of the Wigner function pw and its classical counterpart pc are negligible. 

Independent of the value of Coj we see similar dynamical behavior with identical stages: We start with (i) a rapid 
increase in G, which can be understood as essentially the behavior of a closed quantum system since the gradients 
of the distribution have not increased enough for the diffusive term to become relevant. This is followed by (ii) a 
turnover and an exponential decrease in G as the distribution starts filling the phase-space and the diffusive terms 
kick in - the overall distribution continues to evolve, but the relative size of the quantal terms is now decreasing as 
a result. Finally, (iii) at long times G saturates; this happens because the distribution has almost relaxed to its final 
state. At this point final stage, both Tr[Lq(t)^] and Tr[pu,(t)^] decrease exponentially with time. That is, at longer 
times, the quantal and classical contributions to the evolution of the Wigner function reach a steady-state ratio. 

To examine the functional dependence of G(t) on Coj we start at Co 0, the classical limit for open quantum 
systems. Fig. [2] shows that in the classical limit the greatest value of G(t) for small Co is approximately 

max G(t) w 25 (1 - exp(-0.08Co)) ■ (3) 
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FIG. 2: This plot shows in when we re-scale G{t) by dividing the data for < 2 in Figure[T]by the maximum value of G{t), the 
data collapses to approximately a single curve. Also this data shows that the majcimum value of G{t) as a function is described 
very well by a Gaussian in (o. 



That is, the maximum value of G{t) scales approximately quadratically with in the limit of small Co and is only 
a function of Co- This is a remarkably time-independent relationship, and in Fig. ([2]) we show that for Co < 2, G{t) 
collapses to essentially a single function of Co when we divide through by the maximum value given in Eq. 

The transition out of the near-classical regime starts at Co ~ 2. As C increases further, we expect quantum effects 
to increase and might naively predict that Tr « TT{dtPwit))^ before exceeding it and then becoming the dominant 
term. However, the saturation effects in Fig. ((!]) indicate otherwise. Further, plots (Fig. ^) of the absolute value of 
the Wigner function for the Duffing Oscillator evaluated at the times t = 5 and 20 for Co = 0.2 and Co = 10, clearly 
show that as quantum interference becomes more substantial, the amplitude of interference fringes approaches that of 
the classical phase space structure underlying it. However, these classical structures (in this case the noise-broadened 
homoclinic tangle of the stable and unstable manifolds) always retain a substantial contribution, again pointing to 
the wisdom of considering the classical dynamics when thinking about quantum chaotic systems. 

This saturation effect is examined in greater detail in Fig. [JJ in which it is shown that the observed saturation effect 
occurs because the magnitude of TrL^ becomes comparable to Tri^, while the diffusion term remains negligible. This 
saturation can be understood in the asymptotic regime by considering the influence of quantum interference effects 
on pu,. 

Quantum mechanical effects often appear in physical systems as interference patterns with characteristic wavenum- 
bers that scale with h^^. Therefore we expect that if /D is large we expect that dpPw ~ and dxPw ~ dpPw ~ 
,h^dpPw ^ and Dd^pw ~ ^ for a chaotic system. Hence the quantal terms should be comparable to the 
classical in the limit of large Co if Pw is dominated by interference effects, which is reasonable for our Co = 10 data 
according to Fig. ([3]). 

To summarize, we considered the normalized size of the quantum terms in the Master equation for Wigner function 
evolution in an open quantum system. For a driven Duffing oscillator, this measure shows remarkably precise scaling 
over long time-scales with the parameter Co = Tt^ /D. We also see that, independent of Co the dynamics follows a 
similar pattern. For small Co all of our curves collapses to essentially a single curve when scaled by the maximum 
value of the quantum-classical difference. In both limits of large and small Co we see a saturation effect in the size 
of the quantum-classical difference; that is, the instantaneous difference between quantum and classical evolutions 
cannot be either too small or too large. This further confirms the growing intuition that decoherence softens the 
quantum-classical transition for nonlinear systems. Open questions include whether this remarkable scaling of G{t) 
with Co is exhibited for other Hamiltonians, such as the driven rotor, whose Moyal series does not terminate. 
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FIG. 3: This plot shows the absolute value of p„ for Co = 0.2 {h = 0.01, D = 5x lO"") and Co = W (h = 0.125, D = 1.5625 x 10"^) 
evaluated at times ranging from t — 5 to t — 20. The Wigner function for Co = 10 shows clear signs of quantum interference, 
but the data for Co = 0.2 shows no visible signs of quantum behavior at the times considered. This reflects the observation 
that G{t) is a slowly decreasing function for Co ~ 10 because the quantum effects plotted are also persistent with time. 
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FIG. 4: This plot shows that for Co = 10 with {h = 0.125, D — 1.5625 x 10^"^) the classical terms are roughly comparable to 
the quantal terms. For Co = 0.2 with (h — 0.01, D = 5 x 10"'*) the classical terms are substantially larger than the quantum 
terms. The diffusion term is negligible compared to the classical in both cases. 
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